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VESTIGE OF THE STAR CLUSTER BURST IN M51 
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ABSTRACT 

We present a study of the star cluster formation in M51 based on the image data taken with 
Advanced Camera for Surveys (ACS) and Wide Field Planetary Camera 2 (WFPC2) of Hubble Space 
Telescope (HST) . We have derived the star cluster formation rate using the ages and masses of about 
2,000 star clusters estimated by comparing photometric data (F336JU, F435PL, F555W, and F8UW) 
with theoretical population synthesis models. The star cluster formation rate increased significantly 
during the period of 100 ~ 250 Myr ago. This period roughly coincides with the epoch of dynamical 
encounters of two galaxies, NGC 5194 and NGC 5195, expected by theoretical models. The age 
distribution of the star clusters also shows two peaks at about 100 Myr and 250 Myr ago. The star 
cluster mass ranges from 10 3 to 10 6 A/q and the mass function can be represented by a power law with 
an index ranging from a = -2.23 ± 0.34 for t < 10 Myr to a = -1.37 ± 0.11 for t > 100 Myr. The 
mass function of star clusters older than 10 Myr also appears to display the steepest distribution with 
a w —1.50 at around 200 Myr ago, near the expected epoch of the galaxy interaction. We also confirm 
the correlations of cluster size increasing with cluster mass (with a best fit slope of 0.16 ± 0.02), and 
with cluster age (0.14 ± 0.03). 

Subject headings: galaxies: individual (M51; NGC 5194; NGC 5195) — galaxies: interaction — galax- 
ies: spiral — galaxies: starburst — galaxies: star clusters 



1. INTRODUCTION 

Star clusters are an excellent tool to probe the forma- 
tion and evolution of galaxies due to their high bright- 
ness enough to be detected even in some nearby galaxies. 
They are considered to be composed of stars that formed 
at the same time. Some star clusters are also known for 
their longevity: globular clusters are as old as the Uni- 
verse itself. Therefore, we can investigate the very an- 
cient star cluster formation events as well as the recent 
ones in a galaxy by estimating the age of star clusters 
and constructing their age distribution. 

M51 is an interacting galaxy system that is com- 
posed of a Sbc galaxy (NGC 5194) and a fainter SB0 
galaxy (NGC 5195). M51 is an optimal and interesting 
target for the study of star cluster formation for sev- 
eral reasons. (1) It is abundant with interstellar me- 
dia (see a recent study by Schuster et al. 2007 and 
references therein) that serve as basic ingredients for 
stars and star clusters. HII region studies (e.g., Scov- 
ille et al. 2001, Lee, J. H. 2009 in preparation) showed 
that current formation of stars and star clusters is ac- 
tive ly going on in M51. (2) There are some theoreti- 
cal (|Toomre fc Toomre Hl972tlSalo & Laurikainenl^OOOf ) 
and observational (jDurrell et al. Il2003f) studies suggest- 
ing that the two galaxies in M51 experienced a single or 
multiple encounters some hundred Myr ago. Therefore 
it is possible to find any correlations between these dy- 

* Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained from the data archive at the Space Tele- 
scope Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under NASA con- 
tract NAS 5-26555. 

3 Japan Society of the Promotion of Science Research Fellow 



namical events and the star formation e vents in M51. (3) 
M51 is located at a distance of 8.4 Mpc (jFeldmeier et al. 1 
119971 ) so that, with the resolving power of the Hubble 
Space Telescope (HST), it is possible to resolve, at least, 
some star clusters in M51. 

Several st udies used the HS T data to survey star clus- 
ters in M51. iBik et al. I (|2003h detected about 1000 star 
clusters in a small region near the center of M51 and 
found that the cluster formation rate decreases with age 
during the period of 10 — 1000 Myr ago. However, they 
found no evidence for an increased cluster formation rate 
around 200 — 400 Myr ago, when the fir st encounter of 
NGC 5194 and NGC 5195 was expected. IBastian et al~l 
(2005) showed that the cluster formation rate increased 
about 50 — 70 Myr ago when the second (or last) pas- 
sage of NGC 5195 was expected by a multiple passage 
model ()Salo fc Laurikainen I l2Q00f ) . Still, no clear sign 
of increased star cluste r formation around th e 200 — 400 
Myr ago was shown by IBastian et al. I (|2005D . A hint on 
the increased cluster formation rate around this epoch 
was given by ILee et al. I ()2005D . They derived the age 
distribution of about 400 resolved clusters in M51 using 
the HST archival data and showed that, when compared 
with another typical late type galaxy M101, there are an 
increased number of star clusters with ages of 100 — 500 
Myr, which are consistent with the expected epoch of the 
first encounter. However, their study did not cover the 
entire area of M51 and the number of star clusters used 
in their study was small. 

One limiting factor that has haunted previous stud- 
ies on star clusters in M51 is the limited field of view 
and the insufficient depth of the available data. It be- 
came possible to reduce this problem, when the deep and 
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wide field image data taken with the Advanced Cam- 
era for Surveys (ACS) onboard the HST were released 
(|Mutchler et al. 1120051 ) . Now it enables to investigate the 
star clusters concentrated not only in the disk of NGC 
5194 but also those around the outer halo including the 
companion galaxy NGC 5195 in a homogeneous quality. 
This has led to several new stu dies on the star clusters 
in M51: IHwang fc Lee I (120061) on faint fuzzy clusters 
in NGC 5195. IScheepmaker et al~1 (120071) on the clus- 
ter size distribution. IHwang fc LeeT (l2008|) on the ph oto- 
metric properties of star clusters, lHaas et al. I (120081) on 
the clu ster luminosity function, and IScheepmaker et al. 1 
(2009) on the distribution of star cluster formation in the 
disk. However, there is not yet any detailed study using 
this new data on the correlation between the star cluster 
formation and dynamical interactions in M51, and it is 
not yet clear how the star cluster formation was affected 
by the dynamical events in M51. 

In this study, we investigate the age distribution of 
M51 star clusters and derive the cluster formation rate 
using the cata l og of the star clusters in M51 given by 
IHwang fc Lee I (|2008f l. This star cluster catalog in- 
cludes about 3,600 clusters with Vf555W < 23 mag 
and provides photometric data in F435W, F555W, and 
F8UW bands. We compare the photometric data in 
this catalog with the pop ulation synthesis models by 
iBruzual fc Chariot I (|2003f) to estimate the age and mass 
of star clusters. A preliminary study of the age dis- 
tribution of M51 sta r clusters was already presented in 
IHwang fc Lee I (j2007l ) , implying a possible correlation be- 
tween the dynamical encounters of galaxies and the age 
distribution of star clusters in M51. However, this re- 
sult was based on F435IF, and F8UW band 
photometry data, and could be severely affected by the 
age-reddening degeneracy in age estimation. Therefore, 
we included the F336W band image data of M51 taken 
with HST WFPC2 to improve the accuracy of age esti- 
mation for this study. 

The organization of this paper is as follows: Section 
[2] describes the data. Section [3] briefly introduces the 
star cluster selection and photometry procedures along 
with the reduction of F336W band data. Section [4] de- 
scribes the age estimation method. We present the main 
results including the age and mass distribution, and the 
star cluster formation rate in Section [5l We discuss the 
primary results in Section [6l Finally, a summary and 
conclusion is given in Section [71 

2. DATA 

The main data used in this study are F435W, F555W, 
and F8UW band images taken with HST ACS through 
the Hubble Heritage program 10452 (PI: S. V. W. Beck- 
with) . We also used the F336W band images taken with 
HST WFPC2 for six fields through the observing pro- 
gram 10501 (PI: R. Chandar), and another F336W band 
image set that covered the central region of NGC 5194, 
available in the HST archive. A brief information on the 
^336^ band data set is provided in Table [TJ Figure Q] 
shows the observed fields overlaid on the 12' x 12' Dig- 
itized Sky Survey (DSS) image of M51. Seven WFPC2 
F33QW fields cover a major part of M51, but their cov- 
erage is still much smaller than the ACS field. 

Data reductions for ACS data were carried out by 
the STScI including multi-drizzling and image combi- 
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Fig. 1.— Locations of the fields with the HST ACS data (field 
#8 in dashed line) and WFPC2 data (solid line) used in this study 
overlaid on a 12' X 12' Digitized Sky Survey image of M51. Detailed 
information on each pointing is listed in Table jT] 

nation. More detailed informat ion on these HST AC S 
data reductions was given in iMutchler etaTI (f2005h . 3 
All WFPC2 data that were retrieved from the HST 
archive were reduced using the "on-the-ffy" calibration, 
which automatically uses the best reference files for cal- 
ibration. The WFPC2 data reduction steps include bad 
pixel masking, analog-to-digital correction, bias and dark 
subtraction, and flat-field correction. There are two in- 
dividual images per each WFPC2 pointing. These two 
WFPC2 image data were combined using the combine 
task in STSDAS package of IRAF 4 with CRREJECT 
option in order to eliminate cosmic rays and increase the 
signal-to-noise ratio. 

We adopted a distance of 8.4 ± 0.6 Mpc [(to — M) = 
29.62] to M51 determined from the planetary neb ula lu- 
minosity function in M51 (Fcldmcicr et al. 1 11997). The 
corresponding linear scale is 40.7 pc arcsec -1 . The fore- 
ground reddening toward M51 is low, E{B — V) = 0.035, 
and the corresponding extinctions are Ab = 0.150, Ay = 
0.115 , Ai = 0.067, and A v = 0.190 mag (jSchlegel et al. I 
1998). The total magnitudes and colors are B T = 8.96 ± 
0.06 mag and (B T - V T ) = 0.60±0.01 for NGC 5194, and 
B T = 10.4 5 ±0.07 mag and (B T - V T ) = 0.90 ±0.01 for 
NGC 5195 (|de Vaucouleurs et al. |[T99lh . At the adopted 
distance (without internal extinction correction) , the ab- 
solute magnitudes are Mg = —20.81 mag, My — —21.38 
mag for NGC 5194, and = -19.32 mag and My = 
-20.19 mag for NGC 5195. 

3. STAR CLUSTER SELECTION AND PHOTOMETRY 

We detected and classified star clusters in M51 based 
on HST ACS data and made a catalog of the star clus- 
ters with F555W < 23 mag. Detailed information on 
star cluster detection, photometry, and classification was 

3 See also http://archive.stsci.edu/prepds/m51/ 

4 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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TABLE 1 

A list of HST WFPC2 F336W band data for M51 used in this study 



Field No. Prop. ID 



RA (J2000) 
[hh mm ss] 



Dec (J2000) 
[dd mm ss] 



Instruments Filters 



Exposure Times 
[sec] 



7375 
10501 



13 29 47.81 
13 30 00.52 
13 30 07.80 
13 29 51.68 
13 30 05.90 
13 30 00.87 
13 29 46.32 



+47 11 32.84 
+47 15 40.40 
+47 14 00.50 
+47 15 02.85 
+47 11 40.00 
+47 09 39.37 
+47 08 34.39 



WFPC2 



F336W 
F336W 
F336W 
F336W 
F336W 
F336W 
F336W 



2 x 600 
2 x 1300 
x 1300 
x 1300 
x 1300 
x 1300 
x 1300 



described in lHwang fc Lee 1 (|2008h and only a brief sum- 
mary is provided here. 

We carried out the source detection in the background- 
subtract ed F555W band image of ACS data using SEx- 
tractor (|Bertin fc Arnouti [T996). A detection threshold 
of 4 (j and minimum contiguous detected area of 5 pix- 
els were used for finding sources. The flux of the de- 
tected objects was measured using SExtractor in dual 
mode for the FA35W, F555W, and F8UW band im- 
ages. An aperture with r = 6 pixels (0.3" ~ 12 pc) was 
adopted for the photometry. The instrumental magni- 
tudes in the FA35W, F555W, and F8UW bands were 
calibrated using the photometric zero points of the Vega 
magnitude system for the HST ACS/WFC provided by 
iSirianni et al. I (|2005l ). We visually inspected candidate 
sources with F555W < 23 mag considering their radial 
profile, morphology, and environmental conditions. We 
selected and classified star clusters into two different cat- 
egories: Class 1 star clusters that have a circular shape 
and no prominent nearby neighbors; and Class 2 star 
clusters that have an elongated shape and/or irregular 
structure, and/or multiple neighbors. Finally, we se- 
lected 2,224 Class 1 and 1,388 Class 2 star clusters with 
F555W < 23 mag and their catalog was provided in 
IHwang fc Leel (l200l . 

We carried out an independent source detection and 
photometry in the F336W WFPC2 image data. Sources 
in WFPC2 image data were detected using SExtractor 
with a detection threshold of 3 a and minimum contigu- 
ous detected area of 3 pixels (0".3). The flux of detected 
sources was measured using an aperture with r = 0".3 (6 
pixels for PC and 3 pixels for WF) in SExtractor. The in- 
strumental magnitudes in F336W were calibrated using 
the photometric zero points of the Vega magnitude sys- 
tem for WFPC2 listed in the HST Data Handbook for 
WFPC2. A correction for the charge transfer efficiency 
(C TE) loss was m ade based on the prescription given 
bv iDolphinl (|2000f ). 5 The coordinates of the detected 
sources were derived using metric task in STSDAS pack- 
age of IRAF that corrects for the geometric distortion of 
each chip in WFPC2. The derived coordinates of sources 
in WFPC2 data were calibrated to the system of HST 
ACS data. 

We compared the F336W band catalog with the cata- 
log derived from the F435VF, F555W, and F8UW band 
data. There are about 1,400 Class 1 and about 1,000 
Class 2 star clusters that have counterparts detected 
in the F336W band within a matching radius of 0".4. 



5 See http://purcell.as.arizona.edu/wfpc2_calib/ for further up- 
dated information. 
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Fig. 2. — (B — V)- (U — B) diagrams of the Class 1 star clusters 
(left panel) and Class 2 star clusters (right panel) with V < 23 mag 
in M51. The solid line represents a theoretic al evolutionary track 
with Z = 0.02 by IBruzual fc Chariot I H2003I) shifted according to 
the foreground reddening E{B — V) = 0.035. Numbers that run 
from 6.0 to 10.2 with tick marks on the track represent log(age). 
The arrow in each panel shows the direction of reddening corre- 
sponding to E(B - V) = 0.3. 

We checked the spatial distribution of the sources that 
have no F33QW band counterpart and the locations of 
WFPC2 observation fields to separate star clusters that 
do not emit a significant flux in the F336W band from 
those that were not covered in the F336W band obser- 
vations. It is found that the number of star clusters that 
were not observed in the F336W band is about 300 for 
Class 1 and about 80 for Class 2. Hereafter, we use U, 
B, V, and / to denote F336W, F435W, F555W, and 
F&14W, respectively. 

Figure [2] shows the (B — V) - (U — B) color-color 
diagram for the star clusters (Class 1 and Class 2) in 
M51. We also display a theoretica l evolutionary track 
for th e simple stellar population of Bruzual & Chariot 
(2003) with metallicity Z = 0.02, shifted according to 
the foreground reddening E(B — V) = 0.035. It is seen 
that most star clusters detected in U band are blue with 
(U — B) < 0.2 and younger than 1 Gyr. There are only 
a few Class 2 star clusters with (U — B) > —0.4 (older 
than about 100 Myr in age), while there are still many 
Class 1 star clusters with -0.4 < (U - B) < 0.0. It is 
noted that the (U — B) color for the blue clusters varies 
from -1.8 to -0.2, while the (B — V) color varies only 
from 0.2 to -0.2. This shows the usefulness of U band 
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Fig. 3. — Recovery rates of artificial star clusters as a function 
of magnitudes in V and U band derived from Field A (with higher 
crowding) and B (with lower crowding). Solid lines are polynomial 
fits to the data and dashed lines are straight lines connecting data 
points. 

photometry in estimating the age of young star clusters. 

We tested the completeness of star cluster detection 
in V and U band images. Since the efficiency of star 
cluster detection can be affected by the degree of source 
crowding, we selected two separate fields with different 
source crowding for this test: Field A (a = 13 : 30 : 
03.2, S = +47 : 13 : 32.6) with higher crowding and 
Field B {a = 13 : 30 : 05.8, 6 = +47 : 11 : 49.5) with 
lower crowding. Each field was set to a square with size 
of about 1' x V . The selected two fields were covered 
in both V and U bands. Artificial star clusters were 
generated using mkobjects task of ARTDATA package in 
IRAF. The point spread function (PSF) images required 
for artificial cluster generation were prepared in ACS V 
band data and in WFPC2 U band data. 

We generated 120 artificial clusters in one test field im- 
age over the magnitude range of 18 ~ 23.5 mag with 0.5 
mag step (10 artificial clusters per each magnitude step). 
The FWHM and ellipticity of artificial clusters were de- 
rived from the observed size and ellipticity distr i butio n 
of M51 star clusters in Table 1 of IHwang fe Lee I (|2008f) . 
To do that, we modeled observed distributions using a 
Gaussian centered at 3.25 pixels with a = 2.0 for FWHM 
and another Gaussian centered at 0.15 with a = 0.25 for 
ellipticity. Then, FWHM and ellipticity were randomly 
derived from these Gaussian distributions and assigned 
to each artificial cluster. However, we selected only arti- 
ficial clusters with FWHM > 2.0 pixels to be consistent 
with the observed data. We randomly distributed artifi- 
cial clusters in a test field image and carried out source 
detection and photometry in the same way as it was done 
for M51 star clusters. This routine of test was repeated 
10 times in each test field and the results of source de- 
tection and photometry were compared with the input 
model parameters to calculate the completeness of our 
star cluster photometry data. 

Figure [3] shows the result of the completeness test for 
Field A (high crowding) and Field B (low crowding). It 
can be seen that the recovery rate (completeness) of ar- 



tificial clusters with V < 23 mag is higher than 80% 
(> 84% for Field A and > 93% for Field B). However, 
the recovery rate of artificial clusters with U = 23 mag 
drops to about 40% in Field A and about 80% in Field B. 
This shows that the number of star clusters that we may 
have failed to detect increases more significantly in the 
U band than in the V band as the clusters get fainter. 
However, it is noted that the recovery rate of star clus- 
ters is higher than 80% when star clusters are brighter 
than U = 22 mag (84% for Field A and 96% for Field 
B). This value is consistent with the recovery rate of star 
clusters with V = 23 mag, t he magnitude limit ad opted 
for the star cluster survey in IHwang fc Leel (|200l . 

4. AGE AND MASS ESTIMATION 

We compared the photometric data of M51 star clus- 
ters with the model spectral energy distribution (SED) to 
estimate the ages and masses of the star clusters. Model 
SEDs were deriv ed from the theoretical evolu tionary syn- 
thesis model by iBruzual fc Chariot! (|2003f ). For these 
models we adopted a Salpeter initial mass func tion with 
a power slope of x = —2.35 (|Salpeter I I1955D . a lower 
mass cutoff of 0.1 M & and an upper mass cutoff of 100 
M Q . The models span ages from 1 Myr to 15 Gyr for 
metallicity range of Z = 0.0001 ~ 0.05. However, due 
to the age-metallicity-reddening degeneracy in the inte- 
grated colors of star clusters, we assume the solar metal- 
licity for our working models to derive only age and red- 
dening. Observations of the HII regions in M51 shows 
that the current metallicity of the gas in this galaxy is 
approximately solar (e.g., Diaz et al. 1991; Hill et al. 
1997). It is also known that the influence of metallicity 
becomes more pronounced for old star clusters with age 
of > 1 Gyr in their integrated colors. Because the ma- 
jority of star clusters in M51 is expected to be yo unger 
than 1 Gyr from the study in lBastian et al. I (|2005| ) and 
IHwang fc Lee I ()2008| ). the adoption of the solar metal- 
licity is considered appropriate. 

We compared an SED model with a certain age to 
the observed data for each star cluster, after redden- 
ing the model by E(B — V) values between 0.0 and 
0.6 in steps of 0.015, for about 2,400 star clusters with 
UBVI photometry. This range includes the reddening 
values derived from the previous studies for the most 
star clusters in M51. lLamers et al. I (|2002h showed that 
the reddening in the bulge of NGC 51 94 is typically 
E{B - V) w 0.2 and iLee et al. I (|2005f ) also reported 
that most star clusters in M51 appear to have associated 
reddening E(B — V) < 0.1. For every combination of 
age and E(B — V), we fitted the theoretical SED models 
to the observed photometric data and calculated the \ 2 ■ 
The fit with a minimum \ 2 w &s adopted as the best fit 
for the age and E(B — V) combination. 

The number of finally accepted fits is 1,125 for Class 
1 and 835 for Class 2 star clusters. Then we derived the 
masses for these star clusters by combining the M/Ly 
ratio of the best fit model with the measured V band 
cluster luminosity, extinction, and the distance to M51. 
Figure |4] shows the age versus mass diagram for these 
star clusters. The solid curve indicates the magnitude 
limit of V = 23 mag adopted for star cluster selection in 
IHwang fc~ L cc (2008j). The other lines show the expected 
age and mass of star clusters with different U band mag- 
nitude ranging from 23 to 21.5 mag. This shows that 
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Fig. 4. — Mass vs. age of the star clusters in M51 (circles for Class 
1 and triangles for Class 2). The solid line shows the magnitude 
limit of V = 23 mag used to select star clusters in I H wane Leel 
(2008). Other lines display magnitude limits in the U band. 

star clusters with 1 Gyr are detected only if they are 
more massive than 10 5 2 M Q (brighter than U « 22 mag) 
as indicated in Figure |U Note that the completeness for 
U = 22 mag band is 80% as shown in Figure [3] 

Several features are noted in Figured! First, there is a 
region with few star clusters between logt = 7 and 7.5 (t 
represents an age in units of years). This is considered to 
be an artifact that is caused by the age-dating techniques 
that compare the integrated colors of star clusters with 
the synthesis models. In this age range, the predicted 
colors of star clusters span a very narrow range of colors, 
e.g., S(U - B) w 0.16 mag, 5(B - V) g 0.00 mag, and 
S(V— I) ~ 0.22 mag in the models bv lBruzual fc Chariot I 
(2003). Therefore, even a small photometric error can 
result in a fit to star cluster age younger than log t — 7 or 
older than logt = 7.5. Second, there are two age regions 
of over abundance in star clusters: one with logt < 7 
and the other with 7.5 < logt < 9.0. The over-density 
at logt < 7 is suspected to be, in part, due to the under- 
density at 7.0 < logt < 7.5. Nonetheless, it is still a sign 
that shows the existence of many young star clusters, 
which is consistent with the results of several studies on 
HII regions (e.g., Scoville et al. 2001; Lee, J. H. et al. 
2009 in preparation). The other over-density at 7.5 < 
logt < 9.0 indicates that many star clusters were actually 
formed in this epoch. Third, there is an apparent trend 
that the mass of the most massive clusters increases with 
age. This is known as a 'size-of-sample effect' noted by 
IHunter et al. I (2003). However, it seems that the upper 
limit of this relation is composed of multiple components. 
We discuss this further in Section [(TH 

Similar features in the age and mass distribution were 
also reported in several other studies of M51 star clus- 
ters. For example, the under-density of star clusters in 
7.0 < logt < 7.5 ca n be noted in the res ults given by 
iBastian et al. I (|2005h (their Figure 10) and lGieles et al. I 
d2005|) (theirFigure 1) in combination with the strong 
over-density of star clusters with logt < 7. This is con- 
sistent with the result sh own in Figure [4j However, 
IScheepmaker et aTH (|2009f ) show rather different age and 



Fig. 5. — Age distributions of the star clusters in M51:(a) all 
(Class 1 and 2) clusters, (b) Class 1 clusters, and (c) Class 2 clus- 
ters. Three peaks at logt as 6.7, 8.0, and 8.4 are noted. 

mass distributions in their study of M51 star clusters us- 
ing the same data that a re used in this study. Most star 
clusters in the sample of ISchccpmak er~et al. I (|2009f ) are 
found to be younger than logt = 7.5 and there is no evi- 
dent feature of under-density in 7.0 < lo gt < 7.5, which 
is contrary to the results of other studies (|Bastian et al. I 
120051 : iGieles et all 120051 ) as well as this study. One of 
the possible reasons for this difference may be different 
criteria adopted for star cluster selection, which may lead 
to a different sample of star clusters used for the analysis. 

Another note w orthy point shown in 

IScheepmaker et aTH ( |2009j ) is that different popula- 
tion synthesis models can introduce a systematic change 
to the derived age and mass distributions of star clusters. 
However, the result of their study shown in their Figure 
4 reveals that the age distribution of star clusters with 
logt > 8 is found to be more-or-less consistent with 
one another even in the case that different isochrones 
(Padova or Geneva) and different mctallicities are 
adopted for models. Considering these inconsistencies 
and limits caused by the different analysis methods 
adopted, we focus on the age and mass distributions of 
star clusters with logt > 7.5. 

5. RESULTS 

5.1. Star Cluster Age Distribution 

Figure [5] shows the age distributions of all star clusters 
(including both Class 1 and Class 2), Class 1 star clus- 
ters, and Class 2 star clusters. It is clear that the age 
distributions display two over-dense epochs: one with 
logt < 7 and another with 7.5 < logt < 9.0. The over- 
densities of star clusters in these two periods are evident 
in the age distribution of both Class 1 and Class 2 star 
clusters, as shown in Figures [Hh) and (c). However, 
about 1000 Class 1 star clusters are found in the epoch 
of 7.5 < logt < 9.0, while only about 100 Class 1 star 
clusters are younger than logt = 7. On the other hand, 
for Class 2 star clusters, the number of Class 2 star clus- 
ters younger than logt = 7 is about 230, about one third 
of the number of star clusters with 7.5 < logt < 9.0. 
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Fig. 6. — Age distributions of the star clusters in M51 for different 
U band magnitude limits: (a) U < 23.0 mag, (b) U < 22.5 mag, 
(c) U < 22.0 mag, and (d) U < 21.5 mag. Note that three peaks at 
logt ~ 8.0 and 8.4 are visible clearly in all panels, and the oldest 
peak at logt s=s 8.4 is weak for the brightest limit U < 21.5 mag. 

One interesting point is that there are distinct features 
in the age range of 7.5 < logt < 9.0 in Figure [S£a). The 
most pronounced one is a broad component over logt ks 
7.8 - 8.5 that displays two peaks at log t as 8.0 and logt as 
8.4. The same broad component is also present in the 
distributions of Class 1 and Class 2 star clusters shown 
in Figures E^b) and (c), respectively. However, for Class 

1 star clusters, the peak at logt as 8.4 is more enhanced 
than that at logt as 8.0. On the other hand, the peak at 
logt as 8.0 is relatively more significant in case of Class 

2 star clusters. Other than these, the age distributions 
of Class 1 and Class 2 clusters are consistent with each 
other. Therefore, we consider the age distribution of all 
star clusters including Class 1 and Class 2 clusters for 
the further analysis. 

We investigated the effect of incompleteness to the age 
distribution of star clusters. As shown in Figure 21 the 
incompleteness with different U band magnitudes may 
bias the interpretation of age distribution. Therefore, 
we derived the age distribution of star clusters changing 
the limiting U band magnitude starting from U = 23.0 
to 21.5 mag, and showed the results in Figure [6] Al- 
though the amplitude of histograms in 7.5 < logt < 9.0 
is gradually decreasing as the U band magnitude limit 
gets brighter from U < 23 to U < 22 mag, the overall 
features in the age distribution remain consistent. Even 
when the U band magnitude limit is raised to U = 21.5 
mag, where the completeness is estimated to be about 
90% (see Figure[3]), some hint of peaks at logt as 8.0 and 
logt as 8.4 can be noticed, as shown in Figure ISJd). This 
indicates that the two peaks at logt as 8.0 and logt as 8.4 
in the age distribution of star clusters shown in Figure [5] 
are reliable against any incompleteness problem expected 
from the photometry. 

5.2. Star Cluster Mass Function 

Figure [7] shows the mass functions of all star clusters 
(including both Class 1 and Class 2), Class 1 star clus- 
ters, and Class 2 star clusters in M51. It is found that the 
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Fig. 7. — Mass functions of the star clusters in M51: (a) all (Class 
1 and 2) clusters, (b) Class 1 clusters, and (c) Class 2 clusters (solid 
line histograms). Dotted line and dashed line histograms represent 
very young star clusters with t < 10 Myr, and the intermediate-age 
star clusters with t > 100 Myr, respectively. Solid lines represent 
power law fits (linear fits in the logarithmic scales) for given mass 
ranges. 

mass of star clusters in M51 ranges from 10 3 to 10 6 Af Q 
for both Class 1 and Class 2 star clusters and that the 
overall mass function appears to be separated into two 
components by a dip at « 1O 3 7 M . This dip, separating 
low mass and high mass clusters, is caused by an artifact 
of the SED fitting method over logt = 7.0 - 7.5 as de- 
scribed in Section [4] and the star cluster detection limit 
of V < 23 mag adopted in this study, as shown in Figure 
[H The star clusters with different masses also represent 
star clusters with different ages, as shown in each panel 
of Figure [7J the low mass star clusters are younger than 
10 Myr old (dotted lines) and the high mass star clusters 
are mostly older than 100 Myr old (dashed lines). 

The mass functions of young (< 10 Myr) and 
intermediate-age/old (> 100 Myr) star clusters can 
be represented with a power law over the range of 
3.6 < log(M/M Q ) < 4.1 for young clusters and 5.0 < 
log(M/M©) < 5.7 for intermediate-age/old clusters, re- 
gardless of Class 1 or Class 2. It is noted that the up- 
per limit for intermediate- age/old clusters is marked by 
a weak bump at log(M/M©) ~ 5.8, as shown in each 
panel of Figure [7] We derived a power law index a of 
the cluster mass function as defined in NdM oc M a dM . 
It is found from the sum of Class 1 and Class 2 clus- 
ters that the power law index a for old star clusters 
with 5.0 < Log(M/M ) < 5.7 is -1.37 ± 0.11. For 
star clusters with log(M/M ) < 4.5, where the mass 
function of massive star clusters starts to flatten out, 
the index a, derived using star clusters younger than 
10 Myr, is a = -2.23 ± 0.34. This index a of young 
star clusters is in agreement with the cluster mas s func - 
tion index a = -2.1 ± 0.3 reported bv lBik et al~l (|2003l ) 
but it is slightly low er than a = —1.70 ± 0.08 given by 
iGieles et ahl (|2006D for M51 star clusters. It is, how- 
ever, in agreement with the cluster initial mass func- 
tion (CIMF) index a sa ~2-0 of star clusters in late 
type galaxies (|Zhang fc Fall Ill999t Ide Griis et "aT~ll2003l: 
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The overall shape of the mass function of Class 1 and 
Class 2 star clusters is similar to each other. How- 
ever, the mass function of Class 2 clusters appears to 
be slightly steeper than that of Class 1 clusters. That is, 
the power law index a is —2.13 ± 0.35 for young (< 10 
Myr) and -1.62 ±0.37 for old (> 100 Myr) Class 2 clus- 
ters, while a = -1.03 ± 0.21 for young and -1.27 ± 0.14 
for old Class 1 clusters. Although the indices for old clus- 
ters may be in agreement with each other considering a 
relatively large error associated with Class 2 clusters, the 
difference in the mass function indices for young Class 1 
and Class 2 clusters i s rather significant. As explained in 
iHwang fc Lee I (|2008ft and in Section [3j one big difference 
between Class 1 and Class 2 clusters is their morpholo- 
gies: Class 1 clusters are in circular shape, while Class 
2 clusters have elongated or irregular structures. This 
suggests that many low mass clusters are usually formed 
or embedded in elongated and irregular structures. 

To test the effect of incompleteness on the star cluster 
mass function, we have investigated how different lim- 
iting magnitudes in U and V band affect the overall 
shape and the slope of the mass function of star clus- 
ters. We constructed the cluster mass function using 
star clusters selected with varying limiting magnitude 
conditions in U and V bands ranging from 23 to 21.5 
mag, respectively, and derived a power law index a of 
the mass function over the same mass range as shown 
in Figure The resulting index a for young clusters 
(logt < 7) turns out to be the same under the varying 
limiting magnitudes from J7<23to?7<21.5 mag and 
from V < 23 to V < 22 mag, respectively. Even when 
the constraint is set to V < 21.5 mag, the mass function 
index is a — —1.88 ± 0.57, which is still in agreement 
with a = —2.23 ± 0.34 for young clusters shown in Fig- 
ure [7ja). For old clusters (logt > 8), it is also found 
that the mass function index is a ~ —1.30 ± 0.20 for all 
magnitude constraints in U and V bands. This indicates 
that the power law indices of the young and old cluster 
mass functions derived in this study are affected little by 
the incompleteness problem. This is because the power 
law fits were done using the mass limits set sufficiently 
high enough to avoid severe incompleteness problem. 

5.3. Evolution of the Star Cluster Mass Function 

The variation of the slope of mass function depending 
on star cluster age shown in Figure [7] suggests that the 
index of cluster mass function may be a function of time. 
The cluster mass function for given time is determined by 
the initial cluster mass function, the cluster disruption, 
and the fading below the detection limit, as well as the 
variable cluster formation rate. We investigated the mass 
function of star clusters that formed in different epochs 
and derived its power law index for each epoch using the 
combined set of Class 1 and 2 star clusters, and the result 
is shown in Figure [5] Age ranges were chosen to separate 
star clusters roughly according to different features in the 
age distribution, as shown in Figure [5l 

Figures Oa) through [UJe) show that the upper part 
in the mass function of star clusters in each epoch can 
be represented approximately by a single power law and 
that the slopes of mass function change depending on 
the epochs. The derived power law indices of mass 
function are a = —2.23 ± 0.34 for logt < 7.0 (panel 
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Fig. 8. — Mass functions of the star clusters in M51 for different 
age ranges (including both Class 1 and Class 2): (a) logi < 7.0, (b) 
7.5 < logt < 8.0, (c) 8.0 < log* < 8.3, (d) 8.3 < logi < 8.5, and 
(e) 8.5 < logt < 9.0. In each panel, the power law index of mass 
function is derived by fitting data points marked by circles and the 
solid line represents the best fit. Solid lines represent the power 
law fits to the data (marked by circles) for given mass ranges, and 
the resulting power law index for each age range is given in each 
panel (including both Class 1 and Class 2). 

a) , a = -1.22 ± 0.07 for 7.5 < logi < 8.0 (panel 

b) , a = -1.68 ± 0.33 for 8.0 < logi < 8.3 (panel 

c) , a = -1.81 ± 0.27 for 8.3 < logi < 8.5 (panel d), 
and a = -1.37 ± 0.36 for 8.5 < logt < 9.0 (panel e). 
These results are consistent with those in Figure [7] in 
the sense that the mass function of star clusters younger 
than logt = 7 is steeper than that of older star clusters. 
However, one interesting point in Figure M is that the 
mass function of star clusters with 7.5 < logt < 8.5 gets 
steeper as ages increase. Especially, for star clusters with 
8.3 < logt < 8.5 shown in FigurelUJd), the power law in- 
dex a = — 1.81±0.27 is roughly in agreement with that of 
young star clusters shown in Figure [5fa). It is also noted 
that the relatively shallow mass function of star clusters 
with logt > 8.5 (a = — 1.37 ± 0.36) is consistent with 
that of star clusters with logt > 8.0 (a = -1.37 ± 0.11) 
shown in Figure [7] 

Since the cluster mass function index changes system- 
atically depending on the epoch of star cluster formation, 
we investigated the power law index of cluster mass func- 
tion using star clusters with logt = 7.5 - 8.6 when the 
condition of star cluster formation is expected to change 
drastically. The width of the logarithmic age bin was 
set to 5 logt = 0.20 and 0.15 to keep the number of star 
clusters more than 100 in each logarithmic bin, whenever 
it is possible. The maximum number of clusters in a sin- 
gle bin is about 420 and the minimum is about 70. We 
verified that the resulting cluster mass function in each 
logarithmic age bin can be always represented by a power 
law function and derived the index of mass function. 

Figure [5] displays the power law indices of mass func- 
tions derived from logarithmic age bins with S logt = 
0.20 (squares) and 0.15 (triangles). If the power law in- 
dex was derived using more than 100 star clusters in 
each age bin, then the corresponding symbol is plotted 
in solid line. Otherwise, the symbol is marked in dotted 
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Fig. 9. — Mass function power law index vs. age of the star clus- 
ters in M51 (including both Class 1 and Class 2). The mass func- 
tion index was derived for two logarithmic age bins, 5 logt = 0.15 
(triangles) and 0.20 (squares). Note that the mass function index 
shows a minimum value (the steepest) at logt si 8.3. Approximate 
linear fits are drawn for (t) < 8.3 (dashed line) and logt > 8.3 
(dotted line). 

line. Figure [5] shows that, as we move to the higher age 
domain starting from logt w 7.5, the star cluster mass 
function gets steeper, reaching a « —1.50 at logt « 8.3, 
as shown in the dashed line derived from a simple linear 
fit. Then, the trend is reversed and the mass function 
gets shallower again for logt > 8.4. However, the num- 
ber of star clusters used to derive the index at logt « 8.5 
is smaller than 100 (about 70) and the scatter between 
indices marked by a square and a triangle in dotted lines 
is large. This suggests significant uncertainties due to 
the limited number of massive star clusters available in 
this age domain (see Figure HJ . 

Interestingly, the epoch of logt « 8.3, when the cluster 
mass function index seems to attain its minimum value, 
almost coincides with the age distribution peak at log t ks 
8.4 as shown in Figure [5] This period is also coincident 
with the epoch of dynamical encounter of the two galax- 
ies NGC 5194 and NGC 5195 ex p ected by theoretical 
mode ls (jToomre fc Toomrel Il972t iSalo fc Laurikaineril 
2000). Therefore, the minimum value of cluster mass 
function index around this period may have some impli- 
cations related with the increased number of star clusters 
and the dynamical interactions of the two galaxies in the 
M51 system. However, if the cluster mass function was 
affected and the number of star clusters increased dur- 
ing the course of dynamical interactions of galaxies, then 
there should be some impacts on the star cluster forma- 
tion rate, which we focus on in the following section. 

5.4. Star Cluster Formation Rate 

The cluster formation rate is a good indicator of the 
star formation activity. However, one of the difficulties 
of investigating the star cluster formation history is that 
more star clusters disappear due to the fading and/or 
the disruption as well as the observation incompleteness 
as we analyze older star cluster populations. Another 
difficulty is that, as shown in Figures [8] and [9] the star 



cluster mass function does not appear to be constant 
but to evolve depending on the formation epochs. This 
suggests that the number of faint star clusters that we 
failed to detect may be different in different age bins, 
which makes the estimation of the number of the lost 
star clusters even more complicated. 

Although the slope of the star cluster mass function 
depends on the epoch, the mass function itself can al- 
ways be represented by a power law function, as shown 
in Figured] Therefore, based on the assumption that ev- 
ery star cluster mass function has a power law form, we 
constructed a 'simple model' mass function of star clus- 
ters by extrapolating the power law fitted mass function 
over logt ~ 7.5 - 8.6 period. This extrapolation method 
allows to compensate the number of star clusters lost 
due to the detection limit and the incompleteness prob- 
lem of the observed data. However, this method does not 
correct for star cluster disruption. For the model mass 
function, we set the lower mass limit to 10 3 M and 
the upper mass limit to 10 6 5 Mq, since there is no star 
cluster with more massive than 10 6 5 M in our data, as 
shown in Figure [7] 

Using these simple model and observed mass functions, 
we calculated the number of star clusters per each Myr, 
i.e., the star cluster formation rate, in each logarithmic 
age bin over 7.5 < logt < 8.6. For this investigation, 
two mass limited samples of star clusters were chosen 
by selecting star clusters with log (M/M ) > 5.0 and 
4.6 < log(M/M ) < 5.0. These criteria were defined to 
select star clusters located above the detection limit of 
U = 22 mag (see Figure^]) and the completeness of about 
80% (see Figure E) at logt < 8.8 for log(M/M ) > 5.0 
and at logt < 8.5 for 4.6 < log(M/M ) < 5.0. The 
definition of mass limited cluster samples also helps to 
avoid any possible artificial contribution from the low 
mass part of the model mass function, and enables to 
test whether any discrepancy in cluster formation rate 
exists between different mass ranges. 

Figure [10] displays the age distribution and clus- 
ter formation rate derived from star clusters with 
log(M/M s ) > 5.0 (panels a and c) and 4.6 < 
log(M/M Q ) < 5.0 (panels b and d). Solid lines are from 
the observed data, while dashed lines are from the model 
mass function. The observed data of star clusters with 
log(M/M Q ) > 5.0 reveal a peak in the age distribution 
and an abrupt increase of the cluster formation rate at 
logt k 8.4, as shown in Figures [TOT a) and (c). Then, 
the data exhibit a gradual decrease in the age distribu- 
tion and a rough plateau in the cluster formation rate in 
7.8 < logt < 8.4. Very similar features are seen in the 
age distribution and cluster formation rate of star clus- 
ters with 4.6 < log(M/M ) < 5.0 in Figures ITOTb) and 
(d). The peak at logt s=s 8.4 shown in the age distribu- 
tion and the cluster formation rate is consistent with the 
peak in the age distribution found around logt = 8.4 as 
shown in Figure 

The data from model mass function also display a peak 
at logt « 8.4 in the age distribution shown by dashed 
lines in Figures [TOT a) and (b). However, for star clusters 
with log(M/M ) > 5.0, one slightly different feature is 
a bump-like peak at logt ~ 7.8 in the age distribution 
in Figure [TOTa). Although this feature does not appear 
significant, it seems to translate into another peak of the 
cluster formation rate at logt « 7.8 apart from the one 
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Fig. 10. — Age distributions (a and b) and cluster formation rates 
as a function of age (c and d) for two mass ranges: log(M/M@) > 
5.0 (a and c) and 4.6 < log(M/M ) < 5.0 (b and d). Solid line 
histograms represent the observed data and dashed line histograms 
represent the simple model data. Note that there is seen a peak at 
log t rs 8.4 in the age distributions for both mass ranges and in the 
cluster formation rate derived from the model data for the lower 
mass range. 

at logt as 8.4 in Figure flUT c). For star clusters with 
4.6 < log(M/M Q ) < 5.0, no such a bump-like feature 
is seen in the age distribution and no cluster formation 
rate peak exists at logt w 7.8, as shown in Figures [TUlb') 
and (d). The model mass function data reproduce main 
features exhibited by the observed data, including a peak 
at logt w 8.4 and a gradual decrease in the age distri- 
bution, and a plateau in the cluster formation rate in 
7.8 < log* < 8.4. 

Figure [TO] shows that star cluster formation rate was 
raised at logt w 8.4 and was maintained in that height- 
ened state for somewhat extended period until logt » 
7.8. This is true for star clusters with both log(M/M0) > 
5.0 and 4.6 < log(M/M Q ) < 5.0. However, it is also 
probable that massive star clusters with log(M/M0) > 
5.0 have experienced more than one burst of cluster for- 
mation: one at logt as 8.4 and another at logt as 7.8. On 
the other hand, according to the model mass function, 
star clusters with 4.6 < log(M/M Q ) < 5.0 exhibit one 
strong peak of the cluster formation rate at logt as 8.4. 

6. DISCUSSION 

6.1. Star Cluster Formation and Dynamical 
Interactions of Galaxies 

The age distribution shown in Figure [5] displays two 
peaks at logt as 8.0 and as 8.4. Interestingly, the epochs 
of these two peaks coincide with the dynamical interac- 
tion epochs pre dicted by the multiple enc ounter model 
of M51 system (|Salo fc Laurikainen 2000). However, if 
the size of logarithmic age bin is enlarged from the cur- 
rent value S logt = 0.15 to S logt = 0.2, then those two 
age peaks merge into one broad component ranging over 
logt = 8.0 - 8.4. In this case, it is not straightforward 
to confirm the existence of independent peaks in the star 
cluster age distribution and the possible correlation be- 
tween those peaks and the individual dynamical inter- 



actions, even though it is still clear that a large pop- 
ulation of star clusters were formed during this period 
of logt = 8.0 ~ 8.4. The uncertainty is mostly due to 
the typical error in age estimation based on the broad 
band SEP fit techniq ue, about 20% (|Anders et al. 1120041 : 
Ide Cxriis et alHl27)0l . 

Figure [TO] shows that the star cluster formation rate in- 
creased over the same period of logt = 8.0 ~ 8.4. More- 
over, more than one event of significant cluster formation 
rate increase are apparent at logt as 8.4 and logt as 7.8, 
depending on the mass of star clusters. However, the 
cluster formation peak at logt as 7.8 is slightly appar- 
ent only from the model data for massive star clusters 
with Log(M/M ) > 5. It may be possible that the peak 
at logt as 7.8 was made by some stochastic effect and 
that it may not be a real feature because the peak is not 
so strong and involves a relatively large error. However, 
this recent peak in the cluster formation rate is noted to 
be nearly coeval with the peak in the age distribution at 
logt « 8.0, as shown in Figure [5] 

Except for the SED fit error, the star cluster disrup- 
tion and the evolution of luminosity and mass function 
can affect the star cluster age distribution and cluster 
formation rate. Although many uncertainties are re- 
lated with these, it is relatively we ll known from the 
theor etical and observational studies ijFall 1120011 lLarsen I 
2009) that the luminosity function and mass function 
of young star clusters follow, in most cases, power law 
functio ns with a sing l e par ameter a that defines its 
slope. iHwang fc Lee I ()2008[ ) showed that the luminos- 
ity function of M51 star clusters is fitted by a single 
power law with a = —2.59 ± 0.03, which is in good 
agreement with the results in other studies o f M51 star 
clusters, a = -2.5 ± 0.1 bv iGieles et al. I (|2006h and 



a = -2.53 ± 0.06 bv lHaas et al. I (|2008) 



Figure [8] shows that the mass function of star clus- 
ters with different formation epochs can be also described 
with a single power law, although the mass function in- 
dex seems to change depending on the epochs. Even 
when we select star clusters from a single logarithmic 
age bin to construct a mass function, the resulting mass 
function is found to follow a single power law over a pe- 
riod of logt as 7.5 ~ 8.6. However, Figure [HI shows that 
the mass function appears to evolve depending on the 
cluster formation epoch. That is, the mass function gets 
steeper as star clusters get older until logt as 8.3 when 
the index reaches the minimum value a as —1.50, mak- 
ing the steepest mass function. Then, the mass function 
appears to overturn to become a shallower function as 
star clusters get older than logt as 8.3. 

The change in the mass function slope may be ex- 
plained as a result of the evolution of star clusters. Espe- 
cially, the cluster mass function can be shallower due to 
the selective disruption of low mass star clusters. This is 
one of the reasons why the mass function of star clusters 
older than 100 Myr (a = -1.37±0.11) is flatter than that 
of star clusters younger than 10 Myr (a — —2.23 ±0.34) 
as noted in Section 15.21 However, Figure [9] shows that 
the mass function gets steeper reaching the power law 
index a as —1.50 at logt as 8.3. If we assume that the 
star cluster disruption and evolutionary fading should 
apply to all star clusters consistently regardless of their 
formation epochs, the change in the mass function slope 
suggests that the number of star clusters increased sig- 
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nificantly around the epoch of logt = 8.3. 

What our results show is that many star clusters are 
formed around logt = 8.4 (250 Myr ago) and this period 
is roughly in agreement with the first encounter period 
of NGC 5194 and NGC 5195 predicted by theoretical 
mode ls (jToomre fe Toomre I Il972t ISalo fe Laurikainen I 
2000). This may suggest that the dynamical interaction 
between two galaxies in the M51 system induced star 
cluster formation, implying a strong correlation between 
the active star cluster formation and the dynamical in- 
teractions experienced by the ho st galaxies. This i s also 
in agreement with the result of ILee et al. I Q2005) that 
the number of star clusters with 100 - 400 Myr is signifi- 
cantly larger in M51 than in other late type galaxy M101, 
suggesting a correlation between the increased number of 
star clusters and the dynamical interaction(s) at the sim- 
ilar epoch. 

There is a hint, although weak, for another increase of 
cluster formation rates around logt = 7.8, e.g., as shown 
in Figu re [1070). In a study on the cluster formation 
in M51. iBastian et aT~l (|2005l ) reported the existence of 
relatively recent star cluster bursts in M51: one at ps 6 
Myr (logt ps 6.8) and another at ps 60 Myr (log ps 7.8). 
The latter burst coincides with the probable increase of 
star cluster formation at logt ps 7.8. The former star 
cluster burst is reproduced as a peak at logt ps 6.7 in 
the star cluster age distribution shown in Figure [5l No 
report about the incre ase of star clu s ter fo rmation in 
logt > 8.0 was given in IBastian et al. 1 (|2005f ). However, 
it is noted that a weak hint is apparent in their Figure 
12 that implies a slight increase of cluster formation rate 
at logt ps 8.5 for star clusters with log(M/M ) > 4.7. 

6.2. Star Cluster Disruption and Age Distribution 

Among many studies in the literature, ILamers I (|2008f ) 
described very concisely two major models that are used 
for studies of s tar cluster age d istribution, i.e., the 'Bal- 
timore model' dFall et al 112005ft and the 'Utrecht model' 
(|Boutloukos fe Lamers 1 12003ft . The Baltimore model 
is characterized by the mass independent disruption of 
about 90% star clusters in each age dex up to approx- 
imately 1 Gyr. On the other hand, the Utrecht model 
includes the evolutionary fading and the mass dependent 
dynamical cluster dissolution by tidal effects as well as 
the mass independent dissolution of star clusters in their 
very early stages, which they termed as 'infant mortal- 
ity'. These two models predict different outcomes in the 
distribution of maximum cluster mass in each age dex 
(|Gieles fe Bastian1l2008ft and in the clust er age distribu - 
tion in a log(dN/dt) versus log(f) plane (|Lamers 1 ,2008). 
We use our age and mass data of M51 star clusters to 
test these two models. 

One way to derive the cluster disruption rate is to use 
the relation between mass and age of star clusters and a 
'size of sample' effect that more massive star clusters are 
found in older logarithmic age bins, as shown in Figure 
|4j This effect comes from the assumption that when star 
clusters are randomly selected from a physically deter- 
mined mass function, the larger samples are likely to h ave 
the more massive star clusters. IHunter et al. I (|2003ft de- 
scribed how to derive the CIMF index a, independent 
of cluster disruption and evolutionary fading rates from 
the slope defined using the maximum mass and the age 
of star clusters, assuming a constant cluster formation. 
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Fig. 11. — (a) The 1st (squares) and 3rd (triangles) maximum 
mass vs. age of the star clusters in M51 (including both Class 1 
and Class 2). Note that there is a break at logt »s 8.2. Dashed 
lines and solid lines represent linear fits to the data for the 1st and 
3rd maximum masses, respectively, (b) Cluster formation rate vs. 
age of the star clusters in M51 (including both Class 1 and Class 
2) for two mass ranges: \og(M / Mq) > 5) (solid line histogram) 
and \og(M / Mq) > 4.5) (dotted line histogram). Linear fits to flat 
and steep parts are plotted in dashed and solid lines, respectively. 

iGieles fe Bastiaiil (|2008l ) expanded this by showing that 
this slope is also an indicator for the mass independent 
disruption rate of star clusters when we have information 
about the CIMF. The Baltimore model expects about 
90% star clusters in each logarithmic age bin to disrupt 
regardless of their ma ss, which would result in a flat dis- 
tribution as shown in IGieles fe Bastian I (|2008ft . 

Figure lH shows that, although it is hard to define due 
to large scatters, the upper limit of mass and age distri- 
bution of M51 star clusters appears to be composed of 
roughly two linear components: 7.0 < logt < 8.2 and 
and logt > 8.2. Figure UTTa) displays this distribution 
of the maximum mass (open square) and the 3rd maxi- 
mum mass (filled triangle) of M51 star clusters in each 
logarithmic age bin with a step of <51ogt = 0.25. A linear 
square fit to logt < 8.2 returns a slope of 1.31 ± 0.32 for 
the maximum mass and 1.04 ±0.24 for the 3rd maximum 
mass, while another fit to logt > 8.2 returns — 0.32±1.22 
for the maximum mass and —0.21 ±0.26 for the 3rd max- 
imum mass. If we fit this distribution with a single line 
over the range of 7.0 < logt < 9.5, the slope turns out 
to be 0.72 ±0.23. This value is stee per than 0.4 reporte d 
bv lLee et al l (l^005ft. 0.26 ±0.09 bvlGieles et al. I IpoM ). 
and 0.23 bv IGieles fe Bastian I (|2008ft for M51 star clus- 
ters. 

This shows that the maximum mass distribution for 
young star clusters with log t < 8.2 is steep, while that for 
old star clusters with logt > 8.2 is more or less flat. The 
slope of the steep part is even slightly ste eper than 0.74 
for LM C and 0.69 for SMC reported by IHunter et al~l 
(2003), although the slope of about 0.72 derived over 
7.0 < logt < 9.5 is in good agreement with these. How- 
ever, the steep slope derived from star clusters with 
logt < 8.2 is consistent with those derived from LMC 
(0.96) and SMC (0.93) clusters with Log(t) < 8.0 by 
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iGieles fc Bastian I (120081). On the other hand, for the 



Antennae galaxy. iGieles &: Basti an ( 2008;) show that the 
slope is 0.01 for clusters with Log(i) < 8.0, which is con- 
sistent with the Baltimore model. These results suggest 
that the Baltimore model, at least, does not apply to 
M51 star clusters, considering that the flat distribution 
of maximum mass in log t > 8.2 may be du e to the phys- 
ical u pper limit of M51 star cluster mass (|Gieles et al. I 

Hoop). 

Figure [TlT b) displays the age distribution of the mass 
limited star clusters in a logdN/dt versus logi plane. 
It clearly shows a bend at \ogt rj 8.4 and the exis- 
tence of a plateau over \ogt as 7.6 - 8.4, which is a 
characteristic of the Utrecht model. On the contrary, 
the Baltimore model expects a c onsist ently decreasing 
relation as shown by iFall et al~l ((20051 ) in their Figure 
2. Even when we change the mass limit, the result is 
still consistent for log(M/M©) > 5.0 (solid line) and for 
log(M /Mq) > 4.5 (dotted line). Applying a simple linear 
fit returns a slope of 0.12 ± 0.19 for the plateau (dashed 
line) and —1.95 ± 0.35 for the steep part at log* > 8.4 
(solid line). This slope —1.95 ± 0.35 of the steep part is 
in agreement with the value (ss —1. 7) expected by the 
Utrecht model considerin g the error (lLamers ll2008h. I f 



we apply the analysis of iBoutloukos fc Lamers I (|2003l ) 



relating this slope A to the CIMF slope a, we can derive 
the dynamical disruption parameter 7 from the equa- 
tion A = (1 + a)/j (for a < 0). If we assume the 
CIMF with a = —2.0, then the disruption parameter 
is 7 — 0.53 ± 0.10. This is in good agreemen t with the 
result 0.57 ± 0.10 reported by lLamers et al. I (|2005h for 
M51. Therefore, the age distribution of M51 star clus- 
ters can be described better by the Utrecht model than 
by the Baltimore model and the cluster disruption pa- 
rameter derived from our data is in agreement with the 
literature. 

6.3. Correlation of Cluster Size with Age and Mass 

The mass of a star cluster is expected to be propor- 
tional to the cluster size since a star cluster would form 
from a parent molecular cloud once it reaches a critical 
density. If we assume a constant density for a cluster, the 
cluster size would in crease in propo rtion to M 1 / 3 . How- 
ever, it was found bv lLarsen~l (|2004f ) based on the investi- 
gation of young star clusters in 18 nearby spiral galaxies 
that the correlation between the half light radius r e ff and 
the mass (M) of star clusters returned a shallower slope 
than expected fro m a constant-density assumption. For 
M51 star clusters. iBastian et al. I ((20051 ) tried to find any 
correlation between the cluster size and mass but there 
was no apparent cor r elatio n in their data. More recently, 
iScheepmaker et al. I (j2007l ) also investigated the size dis- 
tribution of star clusters using the HST ACS data and 
they also reported that they did not find evidence for any 
direct relation betw een mass and rad ius of the clusters. 
On the other hand, iLee et al~l (2005) showed the corre- 
lation between the size and mass of star clusters with a 
best-fit slop e of 0.14 ± 0.03 , which is in agreement with 
the result of lLarsenl (|2004h within 1.5<r. 

We have investigated the correlation of size with mass 
of star clusters. The size s of star cluste rs were measured 
using ISHAPE package (L arson 1 fl999Th Detailed infor- 
mation on the size measurement and size distribution o f 
M51 star clusters was described in lHwang fc Lee I (2008). 
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Fig. 12. — (a) Effective radius vs. mass, (b) effective radius 
vs. age, and (c) effective radius vs. mass corrected for the derived 
size-mass relation of the star clusters in M51 (including both Class 
1 and Class 2). Solid lines represent the best power law fits (linear 
fits in the logarithmic scales). 

Figure [T27a) shows the distribution of star cluster sizes 
and masses. The power law relation between these two 
parameters is clearly seen in this plot although with large 
scatters. We derived the best linear fit from this distribu- 
tion with a s lope of 0.16 ± 0.02 , which is consistent with 
the result of ILee et al~l (|2005h at the ~ la level. This 
correlation indicates that the cluster mass distribution in 
M51 is shallower than that expected from the constant 
density assumption. 

Figure [T2T b) shows the distribution of star cluster sizes 
and their ages. It is seen that there is a trend of cluster 
size increasing proportional to the cluster age and this 
is another result of the 'size-of-sample' effect discussed 
in Section RT2l The best linear fit to the data, shown in 
a solid line, gives a slope of 0.25 ± 0.03. However, the 
observed size-age relationship may simply reflect the fact 
that the observed mass limit goes up with age, as shown 
in Figure [U Therefore, in order to check whether there 
is an intrinsic trend between star cluster sizes and ages, 
we correct the derived star cluster mass-size relation and 
fit the data again, as shown in Figure [T27c). It results in 
a linear fit with a sl ope of 0.14 ± 0.0 3, which is steeper 
th an 0.06 ± 0.02 by ILee et al. I (|2005l ) and 0.08 ± 0.03 
by IScheepmaker et al. I ( 20071) . The relation is shallower 
than found before correcting the mass-size relation but it 
still shows that the size of star clusters inc reases with age. 
This r esult is in contradiction to that of Bastian ct al. 
(2005) in which a slight trend of the decreasing size with 
increasing age was shown. 

7. SUMMARY AND CONCLUSION 

We present the age and mass distribution of about 
2,000 star clusters in M51 based on the SED fit of 
UBVI photometric data. The star cluster age distri- 
bution displays two peaks at logt « 8.0 (100 Myr) and 
\ogt Ri 8.4 (250 Myr). These peaks roughly coincide 
with the epochs of dynamical encounters of NGC 5194 
and NGC 5195 predicted b y a multiple encounter model 
(|Salo &: La urikainen"1 l2000| h The star cluster formation 
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rate derived from the age distribution suggests that the 
cluster formation rate increased significantly during the 
period of 100 ~ 250 Myr ago, when the dynamical in- 
teraction in M51 is expected by the theoretical mod els 
(jToomre fc Toomre lll97l iSafo fe Laurikainen ll2000h . 

The mass function of star clusters can be represented 
by a single power law, and the index a appears to change 
depending on the formation epoch: a ss —2.23 ±0.34 for 
t < 10 Myr while a « -1.37 ± 0.11 for t > 100 Myr. In- 
terestingly, the mass function index a changes systemat- 
ically depending on the formation epochs of star clusters 
and the star clusters with logt « 8.3 (200 Myr) displays 
the steepest mass function with a w —1.50 among star 
clusters older than 10 Myr. This also implies an increased 
star cluster formation at the corresponding epoch, possi- 
bly incurred by the dynamical interaction of host galax- 
ies. 

The correlation between the maximum mass of star 
clusters and the logarithmic age displays a steep slope 
of about 1.0 for logi < 8.2, while it is more or less flat 
for \ogt > 8.2. The age distribution of mass limited star 
clusters clearly displays the existence of a plateau over 
logi ss 7.6 - 8.4 and a bend at \ogt « 8.4. The steep 
slope of the maximum cluster mass and the existence 
of a bend in the mass limited age distribution indicates 



that the disruption and evolution of M51 star clusters is 
described better with the Utrecht model than with the 
Baltimore model. The cluster disruption parameter 7 
derive d following the analysis of iBoutloukos fc Lamers I 
(|2003|) is 7 = 0.53 ± 0.10 under the assumption of CIMF 
a — —2. 0, which is in good agreement with the result 
given bv lLamers et al. I (|2005l) . 

We found that the size of star clusters are positively 
correlated with the mass of star clusters with a best fit 
slope 0.14 ± 0.03. It shows that the size of star clusters 
increases with age even after the correction of the mass- 
size relation with a slope 0.16 ± 0.02. 
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